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11 Abstract 

12 

13 We compute a series of Jason-2 GPS and SLR/DORIS-hased orbits using ITRF2005 and 

14 the std0905 standards (Lemoine et al. 2010). Our GPS and SLR/DORIS orbit data sets 

15 span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the 

16 Jason-2 orbit frame translational parameters per cycle by the means of a Helmert 

17 transformation between a set of reference orbits and a set of test orbits. We compare the 

18 annual terms of these time-series to the annual terms of two different geocenter motion 

19 models where biases and trends have been removed. Subsequently, we include the 

20 annual terms of the modeled geocenter motion as a degree- 1 loading displacement 

21 correction to the GPS and SLR/DORIS tracking network of the POD process. Although 

22 the annual geocenter motion correction would reflect a stationary signal in time, under 

23 ideal conditions, the whole geocenter motion is a non- stationary process that includes 

24 secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely 

25 tied to the center of mass (CM) of the Earth consistent with our current force 

26 modeling, whereas GSFC’s SLR/DORIS-hased orbits are tied to the origin of ITRF2005, 

27 which is the center of figure (CF) for sub-secular scales. We quantify the GPS and 

28 SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, 

29 which is assimilated into mean sea level (MSL) error, from the omission of the annual 

30 term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, 

31 currently used by the oceanographic community, only the negligence of the annual term 

32 of the geocenter motion correction results in a - 4.67 + 3.40 mm error in the Z- 

33 component of the orbit frame which creates 1.06 + 2.66 mm of systematic error in the 

34 MSL estimates, mainly due to the uneven distribution of the oceans between the North 

35 and South hemisphere. 

36 

37 Key words: Jason-2, geocenter motion, GPS, SLR/DORIS, mean sea level error 

38 

39 1. Introduction 

40 

41 The origin of the International Terrestrial Reference System (ITRS) is defined to be the 

42 center of mass of the Earth system, including oceans, atmosphere and continental water 

43 (McCarthy and Petit 2004). Ideally, the origin of the International Terrestrial Reference 

44 Erame (ITRE), realization of the TTRS, to which the Jason orbits are referenced, should 

45 coincide with the mean center of mass (CM) of the entire Earth system (Blewitt 2003). 

46 Although, the realization of the reference frame, through the geodetic stations, centered 
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47 in the CM, and the separation from physical processes to which the stations are subject, is 

48 a coupled problem. For example, according to the principal of the conservation of 

49 momentum, the CM has to be a kinematic fixed point, invariant to terrestrial dynamic 

50 processes. However, the redistribution of masses in the Earth system causes geocenter 

51 motion and as such seasonal, annual and trend variations between the CM and the center 

52 of figure (geometric center of the outer surface of the solid Earth) of the solid Earth to 

53 which the actual ITRE is referenced for sub-secular time scales (Dong et ah, 2003; 

54 Blewitt et ah, 2001). Metivier et al. (2010) have found that global ice melting on the 

55 Earth can induce long-term displacements of the geocenter particularly along the Z-axis, 

56 toward the North Pole. They have calculated that the geocenter velocity can reach 0.7-0. 8 

57 mm/yr and is today most probably between 0.3 and 0.8 mm/yr. As such, for the purpose 

58 of accurate geodetic observations, having access to a nearly instantaneous geocenter is 

59 extremely important for those missions that can sense geocenter motions to some extent 

60 but are not good enough to measure it well independently (Wu et al., 2012). Eurthermore, 

61 the CM is directly related to satellite orbital motion and so is the most appropriate choice 

62 to model satellite geodetic measurements (Eritsche et al. 2010), such as altimetry. 

63 

64 Erom the above, and given the required sea level infrastructure stability of 0.1 mm/yr 

65 (Cazenave et al. 2009), geocenter motion of the CM with respect to the CE, ideally, 

66 should also be included in the process of precise orbit determination (POD), which is 

67 based on the site crust-fixed coordinates of GPS, SER, and DORIS stations. This 

68 movement can be thought of as a global degree- 1 loading displacement correction to be 

69 applied to the crust-fixed coordinates of the tracking network in order to reference them 

70 to the CM of the whole Earth (Cerri et al. 2010). Hitherto, the lack of a community 

71 consensus on a geocenter model has not allowed the geocenter to be forward modeled as 

72 part of the Jason altimetry orbit standards (Cerri et al. 2010, Zelensky et al. 2010). 

73 Therefore, our motivation for this investigation arises from the fact that the realization of 

74 an orbital frame for altimetry centered in the CM plays a major role in the definition and 

75 calculation of the rates of global MSE rise. 

76 

77 There have been a number of approaches for the determination of the geocenter motion 

78 models such as: (1) by measuring the translation of a tracking network relative to the 

79 center of the geodetic satellite orbits - the “network shift or geometric approach”; (2) by 

80 observing the deformation of the solid Earth due to the surface mass load - the “degree- 1 

81 deformation approach” (Kang et al. 2009, Eavallee et al. 2006, Blewitt et al. 2001). Dong 

82 et al. (2003) suggest the degree- 1 deformation approach produces more stable geocenter 

83 estimates. Geocenter estimates based on satellite laser ranging (SER), have been 

84 previously reported by Chen et al. (1999). When compared to estimates based on 

85 combined atmosphere, ocean and hydrological model outputs, Chen et al. (1999) found 

86 general agreement at the annual period but little correlation at the monthly timescale. 

87 

88 Here, we use two most recent versions of geocenter models: (1) The UT/CSR RE04 

89 monthly geocenter time-series from Cheng et al (2010), which is based on the “network 

90 approach”. Cheng et al (2010) use SER data from five geodetic satellites (EAGEOS 1& 

91 2, Starlette, Ajisai and Stella) to estimate a 5x5 gravity field along with 3 geocenter 

92 parameters (Tx,Ty,Tz). (2) The Swenson et al (2008) time-series based on the “degree- 1 
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93 deformation approach”. The degree- 1 terms of this model are estimated from a 

94 combination of data from the Gravity Recovery and Climate Experiment (GRACE) 

95 satellite mission and the modeled global atmospheric and oceanic effects to the Stokes 

96 coefficients. 

97 

98 The 1st question we try to answer in this investigation is which of the produced Jason-2 

99 orbit frames based on the GPS or the SLR/DORIS data is more closely centered to the 

100 CM. Eor this purpose, we characterize the spurious signals contained in the translational 

101 parameters of the GPS and SLR/DORIS Jason-2 orbit origin. Then we concentrate our 

102 efforts to the centering of the North/South component, hollowing this investigation, the 

103 suitability of the geocenter motion models and their consistency when used by each 

104 geodetic POD approach, GPS or SLR/DORIS, is examined. Then the error in the MSE 

105 from the omission of the geocenter motion is evaluated. 

106 

107 2. Description of the Jason-2 orbit set 

108 

109 Eor the purposes of this investigation we compute the Jason2 precise orbits using a 

110 dynamic and a reduced-dynamic approach (Bertiger et al. 1994). The current GSEC 

111 Jason-2 dynamic and reduced-dynamic orbits have been computed with GEODYN 

112 (Pavlis et al. 2009) processing GPS (Table 1), and combined SER (Pearlman et al. 2002, 

113 Urschl et al. 2007) and DORIS (Tavernier et al. 2006, Willis et al. 2010) data using the 

114 std0905 standards outlined in Table 1 of Zelensky et al. (2010) and Table 7 of Eemoine et 

115 al. (2010). These standards include the GRACE-derived static gravity field EIGEN- 

116 GE04S1 (Eemoine et al. 2007), the GOT4.7 (Goddard Ocean Tide Model) dynamic tide 

117 model (update to Ray, 1999), forward modeling of atmospheric mass flux using ECMWE 

118 pressure data (Klinker et al. 2000), a GRACE derived time varying gravity model 

119 capturing the annual variation (Euthcke et al. 2006), updated ITRE2005 SER and DORIS 

120 station coordinates using EPOD2005 (Ries, 2008, Euceri and Bianco 2007) and 

121 DPOD2005 (Willis et al. 2009) and updated GPS station coordinates and orbits using 

122 IGS05 (IGSMAIE-5447, Eerland and Bourrassa 2006). The GPS std0905 standards are 

123 outlined in Table 1. The GPS constellation orbits are held fixed to the coordinate set 

124 generated from a least squares (ESQ) fit to the IGS05 sp3 orbits. The GPS station 

125 positions are held fixed to their IGS05 coordinates. The entire GPS antenna phase center 

126 and their associated variations map used are compatible with the IGS05 framework. In 

127 the GPS dynamic orbit solution (gpsdyn) once-per-revolution (OPR) along & cross-track 

128 accelerations parameters are included. Eor the SLR/DORIS dynamic orbits 1 OPR 

129 parameter is estimated every 24-h in a 24-h long arc. Eor the gpsdyn orbits OPR 

130 parameters are estimated every 15-h in a 30-h arc with six hrs of overlap between 

131 adjacent arcs (Melachroinos et al. 2011a). Most of the orbit error due to radiation 

132 pressure is characterized by an OPR signal (Zelensky et al. 2010). This signal is largely 

133 removed upon the estimation of empirical OPR acceleration parameters in the orbit 

134 solution (Colombo 1986). However, complex errors in the radiation pressure model 

135 interact with the estimated empirical OPR parameters to create errors largely in the X and 

136 Y, and to a smaller extent, to the Z components of the orbit with a draconitic period of 

137 118-days (Zelensky et al. 2010). Any SRP mis-modeling is expected to have no effect in 

138 the annual variation of the Z-component origin. 
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139 

140 Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 

141 2008) to cycle 74 (July 2010). Next, we compare our gpsdyn to the SLR/DORIS orbits 

142 and the GPS-based reduced-dynamic (gpsred) orbits (Melachroinos et al. 2011b) for 

143 internal validation purposes. Also for reasons of external validation we use a set of 

144 reduced-dynamic orbits from JPL release-1 la standards (jpllla) (Bertiger et al. 2010). A 

145 reduced-dynamic POD solution is based on the denser and geometrically stronger GPS 

146 tracking data (Bertiger et al. 1994, Luthcke et al. 2003). In our gpsred implementation 

147 OPR along & cross-track accelerations are estimated every 30 min with process noise 

148 standard deviation of 1.0 x 10' m/s and an exponential decay function with a correlation 

149 time of Ihr (Table 1). As demonstrated in Melachroinos et al. (2011a) GSFC’s gpsdyn 

150 and gpsred orbits agree to within 1 cm radially with the SLR-DORIS orbits and those 

151 computed from other analysis centers (JPL, ESA and CNES), thus satisfying the accuracy 

152 requirement of 1 cm proposed by the oceanographic community (Cazenave et al. 2009). 

153 Table 2 summarizes the orbit data sets and their associated acronyms used further in this 

154 study. 

155 

156 3. Spurious signals in the Jason-2 orbit origin 

157 

158 Based on a GPS EEO tracking approach (Kang et al. 2009) we extract the Jason-2 orbit 

159 frame translational parameters per cycle by the means of a Helmert transformation 

160 between a set of reference orbits and a set of test orbits. As in Kang et al. (2009), our 

161 GPS EEO tracking system consists of the GPS constellation orbits (fixed), the GPS 

162 ground station network coordinates (fixed) and the GPS onboard Jason-2 receiver in low 

163 Earth orbit. The set of reference orbits is chosen to be the gpsdyn orbits due to the 

164 stronger ties to the force modeling. The dynamic technique provides an orbit mostly 

165 governed by the dynamic modeling while the reduced dynamic technique provides an 

166 orbit mostly tied to the tracking data. The set of test orbits are the jpllla, SLR/DORIS and 

167 gpsred orbits. 

168 

169 To a certain extent, the estimated orbits should follow the TRE origin, as this is defined 

170 by the tracking stations, and for the GPS-based orbits, also by the GPS constellation 

171 coordinates used. This will greatly depend from the techniques used in the POD 

172 processing. In the case of the gpsdyn and SLR/DORIS orbits, the satellite dynamics are 

173 constrained by physical models. The transition from satellite states at different 

174 measurement times to the state at the solution epoch is furnished by integration of the 

175 equations of motion, which are governed mostly by the forces (dynamics) acting on the 

176 satellite over the time of interest. As such the gpsdyn and SLR/DORIS orbit should 

177 supposedly be closer centered to the CM origin defined by the force modeling used in the 

178 POD. In the case of the gpsred orbits, a dynamic and a kinematic tracking technique are 

179 combined for the better elimination of errors related to force modeling. Essentially, the 

180 gpsred orbits should be closer tied to the ITRE origin defined by the geometry of the 

181 denser GPS tracking data. If we suppose that one of the orbits is centered in the CM and 

182 the other in the CP, then the estimated translation parameters could represent in reality a 

183 geocenter motion based on a geometric approach. An apparent advantage of the satellite 

184 tracking approach to interpret the geocenter motion between CM and CP is that they 
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185 determine the absolute loeation of the CM with respect to the Earth’s surface (Wu et al. 

186 2012). For the purpose of clarity, and as Collilieux et al. (2009) mention, in the sub- 

187 secular time scales that we focus on our investigation, we must keep one thing in mind: 

188 we are only able to investigate possible translational variations (what we call “geocenter 

189 motion”) due to the inaccessible constant between CF and CM. 

190 

191 We perform a least squares spectral analysis on the time-series of the estimated 

192 translational parameters. As previously stated, a 118-day signal is dominant in the X and 

193 Y components with the largest amplitudes of 2.8 mm and 2.3 mm respectively (Fig. la 

194 and lb). Especially in the X-component the largest signal comes from the transformation 

195 of the gpsdyn and SLR/DORIS dynamic orbits. The 118-day signal is the precise 

196 draconitic (beta-prime) period for the Jason satellites and this result supports the earlier 

197 discussion about the remaining orbit error due to solar radiation pressure (SRP) mis- 

198 modeling by Cerri et al. (2010) and Zelensky et al. (2010). 

199 

200 Fig. Ic illustrates that in the Z-component the annual signature has the largest amplitude. 

201 Other signals of lower amplitude appear at the 87-days and 112-days period but their 

202 origin remains unclear. These signals are very close to the 4* and 3’^‘^ harmonics (87-days 

203 and 117-days respectively) of the GPS draconitic year of 351 days (Schmid et al. 2007, 

204 Ray et al. 2008). The largest annual signature results from the comparison of the gpsdyn 

205 to the SLR/DORIS orbits. The comparison of the gpsdyn to the jpllla orbits exhibits an 

206 annual signature of smaller, but still, of non-negligible amplitude. This demonstrates the 

207 presence of residual error in the annual frequency either in the gpsdyn, SLR/DORIS or 

208 jpllla orbits. Does this represent an annual motion between the orbit origins or is it only 

209 related to a GPS SRP-induced orbit mis-modelling error, which usually causes orbit 

210 variations at the GPS draconitic year of 351 days? For that purpose we use the gpsred 

211 orbits since those are less sensitive, by definition, to any dynamic mis-modeling errors 

212 (Melachroinos et al. 201 lb). The comparison of the Z-component from the gpsdyn to the 

213 gpsred orbits does not exhibit any significant signatures (magenta line of Fig. Ic), 

214 especially in the annual term but also not in any other draconitic term. The orbit 

215 differences in Figures la-lc do not show any SRP mis-modeling sensitivity in the Z- 

216 component, but rather only sensitivity in the X and Y components of the std0905 gpsdyn 

217 and the SLR/DORIS orbits, at the Jason-2 draconitic period of 118 days. Furthermore, it 

218 seems that the gpsred and gpsdyn orbits both have consistent Z-centering. Thus, the 365- 

219 days signature in the Z-component between the gpsdyn, SLR/DORIS and jpllla orbits, 

220 suggests that it is largely related to an annual motion of the orbit origin. The question that 

221 we now need to answer is which of the orbits is more closely centered to the CM? 

222 

223 4. The north/south centering of the Z-component 

224 

225 In this section we examine the centering of the orbits and turn our attention to the North- 

226 South behavior of the Z-component. As seen previously, the most significant peak of the 

227 Jason-2 translational time-series, exhibits an annual signature in the Z-component. We 

228 incorporated the 3-dimensionnal annual term from the geocenter motion models of 

229 Swenson et al. (2010) and Cheng et al. (2010) inside the Jason-2 POD process (trends 
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230 and biases removed) as a correction to the a-priori position of the tracking stations 

231 according to Dong et al. (2003): 

232 

233 X^^{t) = X^^{t)+X^cF{t) (1) 

234 

235 Where 

236 

237 = + ^^(0 (2) 

i 

238 Here {t)is the observed position under a CF frame (ITRF), X^^and are position, 

239 velocity at reference epoch , defined under the CF frame, Xf^(t) are various nonlinear 

240 time-dependent deformations due to solid Earth tide, pole tide, ocean tide, mass loading 

241 from atmosphere (not applied here), non-tidal oceans (not applied here), surface ground- 

242 water (not applied here), and other local effects, X^p (t)m this case is the time-dependent 

243 degree- 1 loading displacement correction from the annual fitted term of the two 

244 geocenter models where any offsets and trends have been removed. 

245 

246 For the Z-component, a shift in the station coordinates from the geocenter correction does 

247 not correspond exactly in an equivalent shift in the orbit frame due to the orbit 

248 inclination. As such, if there is a shift dz introduced to all a priori coordinates of the 

249 ground stations due to the annual term of the geocenter motion correction model, the Z- 

250 component Tz of the orbit will be affected by (Morel and Willis, 2005): 

251 

252 Tz = kydz (3) 

253 

254 Where A: is a linear transfer function. 

255 

256 At first, we apply the geocenter motion correction models only in the SLR/DORIS 

257 tracking stations of the POD process according to equations (1) and (2). We then 

258 compute the Helmert translational parameters between the gpsdyn and the SLR/DORIS 

259 orbits. Next, an annual curve to each Z-component time-series is fitted and compared to 

260 the annual fitted curves from Swenson et al. (2010) and Cheng et al. (2010). In the 

261 comparison we also include the jpllla and the gpsred orbits. Fig. 2a illustrates the annual 

262 term of the Z-component (black line) by Swenson et al. (2010) and compares it to the Z- 

263 components of the Helmert transformations between the reference orbit set {gpsdyn) and 

264 the test orbits. The SLR/DORIS Z-component annual terms are plotted with and without 

265 the annual term of the geocenter motion correction (blue and red). Fig. 2b compares only 

266 the annual terms of the SLR/DORIS Z-component to the annual term by Cheng et al. 

267 (2010). The Z-component annual amplitudes from the two models and each set of the 

268 Jason-2 SLR/DORIS Helmert transformations are summarized in Table 3. 

269 

270 From the comparison of the SLR/DORIS orbits without the geocenter motion correction 

271 to the gpsdyn orbits (blue line in Fig. 2a and 2b), we find that the amplitude of the Z- 

272 component annual signature is 2.82 mm. After the introduction of the annual geocenter 

273 motion correction from Swenson et al. (2010) and Cheng et al. (2010) in the POD process 
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274 of the SLR/DORIS orbits, the translational variations of the origin drop down in 

275 amplitude by 2.10 mm and 1.17 mm (red line in Fig. 2a and 2b) respectively. The 

276 reduction from the initial annual signature of 2.82 mm between the SLR/DORIS orbits 

277 and the gpsdyn orbits (Table 3), resulted by both geocenter models, represents 25 % of 

278 the Swenson et al. (2010) and 58 % of the Cheng et al. (2010) annual term. In this sense 

279 the Cheng et al. (2010) model performs the best in reducing the SLR/DORlS-gpsdyn 

280 origin translational variations. The total signal reduced with respect to the modeled 

281 annual geocenter motion correction is 39 % for both cases. As expected, both models 

282 propagate consistently as a correction in the POD process. These results suggest that the 

283 SLR/DORIS orbits are not centered in the CM, whereas, the gpsdyn orbits closely follow 

284 the CM origin consistent with the conservative force modeling as this is realized through 

285 the GPS POD processing of Jason-2. Furthermore, the 7-parameter transformation 

286 between the gpsred and gpsdyn orbits, demonstrates that both orbits sets have a very 

287 consistent Z-origin (magenta line in Fig. 2a). Someone would expect that the gpsred 

288 orbits would follow the CF as this is defined by the geometry of the denser GPS tracking 

289 data, which dominate the reduced dynamic technique. On the contrary the gpsred orbits 

290 do not demonstrate any significant Z-origin motion with respect to the gpsdyn orbits, 

291 which further supports the argument that both orbits are centered closer to the CM. The 

292 annual signature of the transformation between jpllla and GSFC’s gpsdyn orbits, (green 

293 line in Fig. 2a) exhibits amplitude of 1.66 mm. In section 3, we have shown the observed 

294 SRP error does not contribute to the annual signature between the jpllla and gpsdyn. 

295 This fact leaves an open question for further investigation to whether the annual signal 

296 seen in the Z-component of 1.66 mm between the two Jason-2 GPS orbit sets from two 

297 different analysis centers is due to an inconsistency in the conservative force modeling of 

298 the two solutions (time-variable gravity field) or whether is due to an origin motion. 

299 

300 5. Geocenter motion and mean sea level 

301 

302 We have found that the orbit set more closely centered to the CM is the Jason-2 gpsdyn 

303 and gpsred orbits. Thus we have succeeded in answering the question opened in section 

304 3. Despite their dynamic definition, the SLR/DORIS orbits are centered in the CF defined 

305 by the SLR/DORIS network. It’s worth noticing that the SLR/DORIS std0905 orbits used 

306 in this study is the official product currently released and used by the community for 

307 MSL estimations. As such, in this section we will characterize the errors in the MSL 

308 studies that would not be seen by the users of satellite altimetry data, when using GSFC’s 

309 gpsdyn or SLR/DORIS orbits (based on the std0905 standards), without a priori 

310 knowledge of the annual geocenter motion. 

311 

312 Primarily, we must be aware that errors in the southern hemisphere MSL estimates 

313 (which are differences between the satellite altitude, as this is provided by the estimated 

314 orbits, and the radar altimetry data) will have a larger effect due to their statistical over- 

315 representation in the radar observations (Morel and Willis 2005). As such, we expect that 

316 the errors due to the omitted annual geocenter motion over the southern oceans, will have 

317 a greater weight. 

318 
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319 In order to characterize these errors and their propagation over the oceans, we perform 

320 the Helmert tranformation between the orbits of the same technique, GPS or SLR/DORIS, 

321 where in one of the solutions the annual term of the geocenter motion has been applied in 

322 the POD process. 

323 

324 The amplitudes of the propagated signals in the Z-component are illustrated in Table 4. 

325 The corrected for the geocenter motion SLR/DORIS orbits exhibit a noticeable annual 

326 effect in the Z-component of 74 % and 81 % compared to each geocenter model. In the 

327 gpsdyn orbits the annual geocenter correction propagates with a ratio of 16 % and 19 % 

328 with respect to each model (Table 4). The resulted amplitudes are small. As the transfer 

329 function of the origin error is different depending on the technique and processing 

330 scheme, the difference between orbits will include a part, which is proportional to the real 

331 geocenter motion. In the case of the gpsdyn orbits the resulted ratios are very small. 

332 

333 Looking at the geographical distribution of the amplitudes of the radial orbit differences 

334 in Fig. 3, we can note the asymmetry in the North-South direction over water when a land 

335 mask is included. The amplitudes of the radial orbit error are significantly larger in the 

336 case of the SLR/DORIS orbits and can reach 2.5 mm in high latitudes depending from the 

337 model. The gpsdyn orbits suffer the least from the omission of the geocenter correction. 

338 As the transfer function of the origin error is different depending on the technique and 

339 processing scheme, the difference between orbits may include a part that is proportional 

340 to the real geocenter motion. In the case of the gpsdyn orbits where the geocenter 

341 correction has been applied, the resulted ratios are small. The above argument supports 

342 the fact that the gpsdyn orbits are indeed centered closer to the CM. Those remain 

343 practically insensitive to the geocenter motion correction introduced as a degree- 1 

344 loading displacement correction to the tracking stations. Which shows that the real 

345 geocenter motion left in the GPS Jason-2 orbits is small since the transfer function 

346 resulted after the transformation is also small (see Table 4). 

347 

348 Fig. 4 illustrates the phases of the geocenter motion correction as this propagates over the 

349 globe into the orbit’s radial component. Even though it is small, it worth noticing that the 

350 GPS technique provides a geographical representation of the propagation of the geocenter 

351 motion into the radial component, similar to the one from Blewitt et al. (2001). 

352 

353 Fig. 5 represents the geographical distribution of the POD omission error on the MSL (in 

354 mm) resulting from the geocenter motion model of Cheng et al. (2010) in the SLR/DORIS 

355 stations for cycle 058 (Jan 28-Feb 07, 2010). The systematic error from the modeled 

356 geocenter motion in the Jason-2 SLR/DORIS orbit frame results in a mean Z-component 

357 of -4.67 + 3.40 mm. This affects the MSL (DH) by 1.06 + 2.66 mm (Table 5). The 

358 systematic error in the Jason-2 gpsdyn orbit frame results in a mean Z-component of only 

359 -0.83 + 0.28 mm which affects the MSL by 0.17 + 0.37 mm. We calculated for both 

360 Jason-2 gpsdyn and SLR/DORIS orbits and compared to previous studies, the functions 

361 (ratios DH/Tz in Table 5) that would result from the Tz error in the orbit frame (and not 

362 in the station’s TRL), in an error in the MSL, due to the negligence of the geocenter 

363 motion. In the case of Beckley et al. (2007), by taking into account k = 0.74 from 

364 equation (4) of Morel and Willis (2005), the real Tz error in the TOPEX orbit frame 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 


365 results into an orbit drift of 1.33 mm/yr resulting from a TRF drift of 1.8 mm/yr in the 

366 stations due to the transition to ITRF2005 from ITRF2000 (Altamimi et al. 2007). This 

367 propagates in a -0.26 + 0.72 mm/yr error in the MSL, which provides a transfer function 

368 of -0.20, following our approach. In Morel and Willis (2005) a Z shift of 10 mm in the 

369 TRF stations is found to propagate linearly in the orbit frame Z-component by 74 % (= 

370 0.74 mm/mm). Given our approach, the error in the Z-component of the orbit frame 

371 resulting from the 10 mm shift would then be 7.4 mm, which results in an error of 1.21 

372 mm in the MSL and a transfer function of -0.16. Our transfer function is closer to 

373 Beckley et al. (2007) and differs by 0.05 only from the transfer function of Morel and 

374 Willis (2005). It seems that all three results are consistent. We should point out that it is 

375 the first time that the study of Morel and Willis (2005) has been verified by real case 

376 scenarios. The tiny differences are probably related to the fact that both in Beckley et al. 

377 (2007) and our investigation the geographical latitudes covered by the real orbit 

378 inclination are not the same with the simulated case scenario analyzed by Morel and 

379 Willis (2005). 

380 Fig. 6 illustrates the observed geographical MSL trend resulting from the geocenter 

381 motion model applied in the SLR/DORIS stations from Cheng et al. (2010) over Jason-2 

382 cycles 001 to 074 (2 years). For the whole period of the SLR/DORIS orbits, the 

383 negligence of just the annual term of the geocenter motion correction creates an apparent 

384 MSL rise of 0.14 + 0.35 mm/yr in 2 years. This is a very important result because if 

385 orbits based on SLR/DORIS and GPS are used during the inter- mission calibration phases 

386 of TOPEX/Poseidon, Jason-1, Jason-2 and future Jason-3 then the omission of the 

387 geocenter motion correction could potentially affect the inter- mission calibration of the 

388 altimeter data . 

389 Finally, by editing the SLR weighted residuals after the application of the geocenter 

390 motion correction we see a slight improvement with the annual term of the Cheng et al. 

391 (2010) model over not using any geocenter motion in the SLR/DORIS orbits. For the 

392 gpsdyn orbits the residuals remain unchanged. 

393 

394 All the above, in addition to the small improvement in the SLR weighted residuals, 

395 demonstrate that the geocenter motion correction should be included as a standard in the 

396 Jason-2 POD process. 

397 6. Conclusions 

398 In conclusion, we have characterized the spurious signals contained in the origin of the 

399 Jason-2 orbits from an analysis of GSFC’s SLR/DORlS-hd&QA and GPS-based dynamic 

400 and a set of reduced-dynamic orbits. A 118-day dominant signal was found in the X and 

401 Y components of the comparison gpsdyn to the SLR/DORIS orbits. The 118-day signal is 

402 the precise draconitic (beta-prime) period for the Jason satellites and it is due to solar 

403 radiation pressure (SRP) mis-modeling. This result is consistent to the analysis by 

404 Zelensky et al. (2010). In the Z-component the annual signature has the largest amplitude. 

405 Other signals of lower amplitude appear at the 87-days and 1 12-days close to the 4th and 

406 3rd harmonics of the GPS draconitic year of 351 days (Schmid et al. 2007), but their 

407 origin still remains unknown. We’ve shown that the comparison of the gpsdyn to the 
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408 SLR/DORIS orbits exhibits a large annual signal in the Z-component suggesting a motion 

409 of the origin between the two orbit sets. The jpllla orbits (Bertiger et al. 2010) also 

410 exhibit an annual signature in Z when compared to the gpsdyn orbits but of smaller 

411 amplitude. This left an open question for further investigation with respect to the 

412 consistency in the force modeling and the origin of the two analysis centers GPS orbit 

413 sets. 

414 

415 We examined the centering of the SLR/DORIS orbits with respect to the gpsdyn orbits 

416 after the introduction of the annual geocenter motion as a degree- 1 loading displacement 

417 correction in the stations. For the geocenter motion correction we have used two models 

418 from Cheng et al. (2010) and Swenson et al. (2010). After the introduction of the annual 

419 geocenter motion correction from Swenson et al. (2010) and Cheng et al. (2010) in the 

420 POD process of the SLR/DORIS orbits, the initial translational variations of the origin 

421 with respect the gpsdyn orbits, dropped down from 2.82 mm in 2.10 mm and 1.17 mm 

422 respectively. Furthermore, the 7-parameter transformation between the gpsred and 

423 gpsdyn orbits, demonstrated that both orbits sets have a very consistent Z-origin. Based 

424 on these facts, we have concluded that our gpsdyn orbits closely follow the CM 

425 consistent with our conservative force modeling, while the SLR/DORIS are centered 

426 closer to the origin of the ITRF, which is the CF for sub-secular scales. Moreover, our 

427 investigation suggests that any SRP mis-modeling error on Jason-2 is not responsible for 

428 the annual signature seen in the Z-component comparison between jpllla and the 

429 GSFC’s gpsdyn orbits. This fact left an open question of whether this annual signal is due 

430 to some difference in the POD modeling, such as the time variable gravity, or whether is 

431 due to an origin motion as is proven to be the case with GSFC’s SLR/DORIS orbits. 

432 

433 We have characterized the errors that would be seen by the users of satellite altimetry 

434 data when using the SLR/DORIS std0905-hased Jason-2 orbits without a priori 

435 knowledge of the geocenter motion. The SLR/DORIS orbits, in which the annuals term of 

436 the geocenter motion has been taken in to account, exhibited a noticeable annual effect in 

437 the Z-component of the orbit frame of 74 % and 81 % compared to each geocenter 

438 model. In the case of the gpsdyn orbits this effect was found to be only 16 % and 19 % 

439 respectively with insignificant amplitude. The gpsdyn orbits remained practically 

440 insensitive to the applied geocenter motion correction. We have depicted that in the case 

441 of the SLR/DORIS orbits the geographical amplitude of the mean radial orbit error (DH) 

442 is significantly larger and can reach 2.5 mm in the poles. Moreover, our transfer function 

443 that connects the error in the Z-component of the Jason-2 orbit frame, from the omission 

444 of the annual geocenter motion correction, to the MSL error of both the gpsdyn and 

445 SLR/DORIS orbits is closer to the transfer function re-calculated from the study of 

446 Beckley et al. (2007) and slightly differs from the one from Morel and Willis (2005). It is 

447 worth noticing that it is the first time where the study of Morel and Willis (2005) has 

448 been revisited with real case scenarios 

449 

450 In this study we addressed only the annual term of the modeled geocenter motion as a 

451 degree- 1 loading displacement to the tracking stations that participate in the POD 

452 process. Indeed, seasonal geocenter motion results from mass transfer at the Earth's 

453 surface. We do model a big part of the degree- 1 signal, but this is only a portion of the 
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454 degree- 1 deformation since we also have associated deformation. The tracking network 

455 displaces also because of higher degrees. Future work could focus on the forward 

456 modeling of the seasonal displacements at the stations together with the complete 

457 geocenter model correction. 

458 

459 We have found that the omission of just the geocenter annual term can contribute to an 

460 apparent 0.14 mm/yr error in the MSL estimates in 2 years based on the std0905 

461 SLR/DORIS orbits. However, although the annual term of the geocenter motion could 

462 reflect a stationary signal in time, the whole geocenter motion is a non-stationary process 

463 that includes secular trends. For example global ice melting on the Earth has been found 

464 to induce long-term displacements of the geocenter particularly along the Z-axis, toward 

465 the North Pole. Metivier et al. (2010) have calculated that the geocenter velocity can 

466 reach 0.7-0. 8 mm/yr and is today most probably between 0.3 and 0.8 mm/yr. Especially 

467 in the last decade it seems that there’s an increase in the geocenter velocity not superior 

468 to 0.5 mm/yr. Since one of the main objectives in the present development of altimetry 

469 MSE is stability at the 0.1 mm/yr level (Cazenave et al. 2009), it would be very 

470 interesting to extend the current study to the whole period of Jason- 1 and Jason-2 with a 

471 complete geocenter motion correction. Our results have shown that the Jason-2 GPS and 

472 SLR/DORIS orbits respond differently to the omission of an annual geocenter model in 

473 the POD process. Hence, if orbits based on SLR/DORIS and GPS are used during the 

474 intermission calibration phases (e.g. TOPEX vs. Jason- 1; Jason- 1 vs. Jason-2), then the 

475 geocenter model omission error could potentially affect the intermission calibration of 

476 altimeter data. We need to elucidate whether this conclusion still applies if other 

477 techniques are used to process GPS data than the ones we have applied in this paper 

478 using fixed and filtered IGS orbits, and GPS double differences on Jason-2. 

479 
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635 Table 1: GSFC Jason-2 GPS POD model standards: std0905 


Reference Frame and displacement of reference points 

GPS 

38 IGS05 (Eerland and Bourassa 2006) TRE stations 

Tidal CoM and 
EOP 

GOT4.7 (update to Ray, 1999); VLBI high frequency 
terms 

Ocean Loading 

GOT4.7 (update to Ray, 1999) all stations 

Earth tide 

IERS2003 

EOP 

lERS Bulletin A daily consistent with ITRE2005 
(Altamimi et al. 2007) 

Precession 

Nutation 

IAU2000 

Satellite surface forces and attitude 

Albedo/IR 

Knocke-Ries-Tapley (1988) 

Atmospheric drag 

MSIS86 (Hedin, 1987) 

Radiation pressure 

8-panel, CR=0.916 (tuned) 

Attitude 

Quaternions 

Tracking data and parameterization 

Tracking data 

Double Difference EC iono-free tracking data, float 
ambiguities, fixed and filtered IGS sp3 orbits 

Troposphere 

1/hr scale(wet-i-dry) troposphere (GME (Boehm et al. 
2006)/GPT (Boehm et al. 2007)-Hopfield) 

modeling 

adjusted using 2 paths (1 station -i- 2 GPS s/c) during the 
POD 


arc : 24-1-6 h long 

Parameterization 

gpsdyn : drag l/8hr, 2 OPR along&cross / arc 


gpsred : drag 1/arc, OPR along&cross/ 15 min, sigma 
=l.e-09, correl time =3600 sec 

Antenna Reference 

GPS stations -i- 

PCOs and PC Vs : igsOS.atx 
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satellites 

JPL GPS antenna PCV map consistent with igsOS.atx, 636 


Revised LC GPS antenna PCO offsets 

GPS Jason-2 

637 

^ 
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642 

643 

644 

645 

646 

647 

648 

649 Table 2: Description of the Jason-2 orbit solutions used in this study 


POD name used in the text Description 

GSFC’s GPS dynamic 

GSFC’s GPS reduced dynamic 

GSFC’s GPS dynamic -i- Cheng et al. 2010 
CoM correction 

GSFC’s GPS dynamic -i- Swenson et al. 2008 
CoM correction 

GSFC’s SLR/DORIS dynamic 

GSFC’s SLR/DORIS dynamic -i- Cheng et al. 
2010 CoM correction 

GSFC’s SLR/DORIS dynamic -i- Swenson et al. 
2008 CoM correction 

JPL’s release- 1 la GPS reduced dynamic 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 


gpsdyn 

gpsred 

gpsdyn_com_csr_ 

gpsdyn_com_swn_ 
SLR/DORIS 
SLR/DORIS _com_csr 

SLR/DORIS _com_swn 
jpllla 
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661 Table 3: Z-component annual amplitudes (mm) from each geocenter motion model and 

662 orbit transformations compared to the ratios of reduction in the annual signature to each 

663 geocenter motion model and the SLR/DORIS - gpsdyn comparison. 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 


Geocenter model 

Annual 

Amplitude 

Geocenter 

model 

applied 

Ratio of the 
reduction to 
the 

SLR/DORIS 
- gpsdyn 
signal 

Ratio of 
the 

reduction 
to each 
model 

Swenson et al. 2010 

1.85 

Cheng et al. 2010 

4.24 

Helmert 

transformation (ref. 
orbit gpsdyn) 

Annual 

Amplitude 

SLR/DORIS 

2.82 

SLR/D ORIS _com_s'^n 

2.10 

Swenson et 
al. 2010 

25 % 

39 % 

SLR/D ORIS _com_c^x 

1.17 

Cheng et al. 
2010 

58 % 

39 % 


664 

665 

666 Table 4: Z-component annual amplitudes (mm) and ratios from each orbit solution after 

667 the geocenter motion correction to each model 


Helmert transformation 

Ref. Orbit 

Amplitude 

(mm) 

Phase 

(degrees) 

Ratio of 
the resulted 
signature 
to each 
model 

gpsdyn_com_swn 

gpsdyn 

0.3 

62.7 

16% 

gpsdyn_com_csr 

gpsdyn 

0.8 

4.9 

19% 

SLR/DORIS_com_csr 

SLR/DORIS 

3.1 

12 

74% 

SLR/DORIS_com_swn 

SLR/DORIS 

1.5 

64.3 

81 % 


668 

669 Table 5: Effect observed on the derived mean sea level (DH) resulting from the Cheng et 

670 al. (2010) geocenter motion correction in the gpsdyn and SLR/DORIS stations for Jason-2 

671 cycle 058 (Jan 28-Feb 07, 2010) 


Orbit comparisons 

Ref. Orbit 

Tz (mm) 

DH (mm) 

DH/Tz 

SLR/DORlS_com_csx 

SLR/DORIS 

-4.67 + 3.40 

1.06 + 2.66 

-0.22 

gpsdyn_com_csr 

gpsdyn 

-0.82 + 0.28 

0.17 + 0.37 

-0.21 
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677 Figure 1: Periodogram (in mm) of the orbit origins after a 7 -parameter Helmert 

678 transformation between the NASA GSFC Jason-2 GPS-based dynamie orbits and the 

679 three test orbits: NASA GSFC Jason-2 GPS-based reduced-dynamie (gpsred), NASA 

680 GSFC Jason-2 SLR/DORIS dynamic and JPL Jason-2 GPS-based reduced-dynamic 

681 (jpllla). a) X-component, b) Y-component, c) Z-component. In purple, blue and orange 

682 are the comparisons to GSFC’s gpsred, SLR/DORIS dynamic orbits and jpllla GPS- 

683 based reduced dynamic orbits respectively. 
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687 Figure 2: Jason-2 fitted annual signal of the Z-component time-series from the 7- 

688 parameter transformation between the gpsdyn (reference orbit) and the test orbits: 

689 SLR/DORIS, jpll la and gpsred. a) compared to the Swenson et al. (2010) applied only in 

690 the SLR/DORIS orbits, b) compared to the Cheng et al. (2010) applied only to the 

691 SLR/DORIS orbits 
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698 


699 





700 Figure 3: Amplitude (in mm) of the geocenter motion correction as it maps into the radial 

701 orbit differences (DH) of the gpsdyn (up) and the SLR/DORIS (bottom) orbit frame. Left 

702 from Cheng et al. (2010) and right from Swenson et al. (2010). 
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Figure 4: Phase (in degrees) of the geocenter motion correction as it maps into the radial 
orbit differences (DH) of the gpsdyn (up) and the SLR/DORIS (bottom) orbit frame. Left 
from Cheng et al. (2010) and right from Swenson et al. (2010). 
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711 Figure 5: Observed geographical MSL error (in mm) resulting from the geocenter motion 

712 model of the SLR/DORIS stations from Cheng et al. (2010) for Jason-2 cycle 058 (Jan 28- 

713 Feb 07, 2010) 
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Figure 6: Observed geographical MSL trend (in mm/yr) resulting from the geocenter 
motion model of the SLR/DORIS stations from Cheng et al. (2010) over Jason-2 cycles 
001 to 074 (2 years) : 0.14 + 0.35 mm/yr 


Detailed Response to Reviewers 


Rev. 1 

While it appears necessary, as demonstrated in this paper, that geocenter motion 
should be applied in the POD process, so that all techniques provide orbits that 
are internally consistent, the authors may wish to consider (or pose to the 
altimetry community), whether the resulting orbits (presumably now properly 
centered at the CM) should in fact then be translated to the CF frame for 
purposes of computing sea level. The question arises since sea level is generally 
measured relative to the crust (e.g., tide gauges). Stated alternatively, should 
orbits used for sea level analysis really be in the CM frame or should they 
translated to the CF frame? This question would apply to the non-tidal as well as 
shorter period tidal geocenter motion. 

Indeed, the reviewer is posing a very interesting and challenging question 
concerning modern and current geodetic sea-level investigations. Our future 
investigation is concentrating at exactly trying to answer to this question with 2 
abstracts and one proposal recently submitted. 

A tide gauge directly measures the displacement of the sea surface relative to a 
land point, and so would seem ideal. Flowever the use of tide-gauge data alone 
to infer global measures of sea-level change is fundamentally problematic due to 
processes that intervene in the relative motions of the sea surface and the solid 
Earth on which the tide gauges are, over a broad range of spatial and temporal 
scales. The land beneath the tide-gauges is subject to motions such GIA (secular 
trends), coastal erosion, sedimentary loading, subsidence, atmospheric loading, 
tectonic processes and the different ocean/land response to present-day mass 
redistributions such as cryospheric loading and terrestrial hydrological loading. 
Another point is that the land motion near the costs on which the tide gauge are, 
provides a very poor sample distribution of the Earth deformation processes that 
will generally not average out on the global scale. Secondly, the vertical land 
motion from geodetic techniques used for the correction of the relative sea-level 
trends from the tide gauge sites, are subject to terrestrial reference frame errors. 
Also, an error in the terrestrial reference frame origin at the center of mass of the 
Earth implies an error in the height of sea surface inferred by satellite altimetry 
observations. Orbits whose origin is closer centered at the center of Mass of the 
Earth system, would ideally be insensitive to those reference frame realization 
errors. Tide gauges are immune to terrestrial reference frame errors only if the 
vertical velocities used for the correction of the relative sea level rates are 
inferred by other technique (cf. internal rates by G. Mitchum studies) and thus 
can be used only locally as the ground truth of calibration/validation for the 
satellite altimetry. In order though for these two types of sea-level observations to 
be comparable they need to be defined in the same ref. frame origin. For the tide 
gauges the origin is the center of the Solid Earth (CE) where as those of the 
satellite altimetry “ideally” is the center of mass of the Earth system (CM) defined 
by the satellite orbit dynamics. However the CE is never a frame of reference that 
can be realized by space geodesy. In GPS practice, for example, the CF frame is 
most commonly used. Also we must not forget that in order for the two types of 



sea-level observations to be comparable, the vertical motions (from GPS heights 
normal to a geocentric ellipsoid) must be removed from the tide gauges records. 
So the way this vertical motion is defined is very important for the inter- 
comparisons, which means improved reference frames that for example would 
take into account seasonal mass redistribution, its effect on degree-1 
deformation and therefore its effect on the frame origin, are needed. 


Rev. 2 


All suggested changes have been applied. 



